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dimensional radial Schrodinger equation with some molecular potentials 
like pseudoharmonic and modified Kratzer potentials are obtained. The re- 
striction on the parameters of the given potential, 5 and -q are also given, 
^3 \ where r] depends on a linear combination of the angular momentum quan- 

tum number £ and the spatial dimensions D and 5 is a parameter in the 
. ansatz to the wave function. On inserting D = 3, we find that the bound 

^ , state eigensolutions recover their standard analytical forms in literature. 
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1 Introduction 

The solution of the fundamental dynamical equations is an interesting phenomenon 
in many fields of physics and chemistry. To obtain the exact Estate solutions of 
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the Schrodinger equation (SE) are possible only for a few potentials and hence ap- 
proximation methods are used to obtain their solutions. According to the Schrodinger 
formulation of quantum mechanics, a total wave function provides implicitly all 
relevant information about the behaviour of a physical system. Hence if it is ex- 
actly solvable for a given potential, the wave function can describe such a system 
completely. Until now, many efforts have been made to solve the stationary SE 
with anharmonic potentials in three dimensions and two dimensions [1-5] with 
many applications to molecular and chemical physics. The study of the SE with 
these potentials provides us with insight into the physical problem under consid- 
eration. However, the study of SE with some of these potentials in the arbitrary 
dimensions D is presented in (cf. Ref.[6] and the references therein). Further- 
more, the study of the bound state processes is also fundamental to understanding 
of molecular spectrum of a diatomic molecule in quantum mechanics [7]. The 
purpose of this paper is to carry out the analytical solutions of the D-dimensional 
radial SE with some diatomic molecular potentials by making a suitable ansatz to 
the wave function. These molecular potentials can be taken as pseudoharmonic 

potential [8,9] V(r) = D e (y- — J , and a modified Morse or Kratzer-Fues po- 
tential [10-15] which is of a Mie-type form (cf. Ref. [14] and the references 

therein) V(r) = D e (j—^-^j where D e is the dissociation energy between two 
atoms in a solid and r e is the equilibrium intermolecular seperation. Morse poten- 
tial [10] is commonly used for anharmonic oscillator. However, its wavefunction 
is not vanishing at the origin. On the other hand, the Mie-type and also the pseudo- 
harmonic potentials [8,9] do vanish. The Mie-type potential [14] has the general 
features of the true interaction energy [16], interatomic and dynamical properties 
in solid-state physics [17]. The pseudoharmonic potential may be used for the 
energy spectrum of linear and non-linear systems [18]. The pseudoharmonic [8,9] 
and Mie-type [13-15] potentials are two exactly solvable potentials other than the 
Coulombic and anharmonic oscillator [9]. This paper is organized as follows. In 
Section [2j we solve analytically the D-dimensional radial Schrodinger equation 
with the molecular pseudoharmonic and modified Kratzer-Fues potentials. We 
also obtain the bound state eigensolutions of these molecular diatomic potentials 
by making a suitable ansatz to every wave function. The results are used to find 
energy spectra of diatomic N 2 , CO, NO and CH molecules [15,19-20]. The 
concluding remarks will be given in Section [3j 
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2 The D-Dimensional Radial Schrodinger Equation 
with Some Molecular Potentials 



In the .D-dimensional Hilbert space, the radial wave function ip(r) of the Schrodinger 
equation can be written as [21] 

d 2 (D-l)d £{£ + D-2) 2/i /rn ir . , n 

d? + -^Tr~ L + $&-V(r))\i,(r)=0, (1) 

where £ denotes the angular momentum quantum number and ji is the reduced 
mass. Furthermore, to remove the first derivative from the above equation, we 
define a new radial wave function R{r) by means of equation 

ij)(r) =r- {D -^/ 2 R(r), (2) 

which will give the radial wave function R(r) satisfying the wave equation 

{£ _ mi+D-wD-mD-w + p E _ v(r)) j m = . 

(3) 

Further, Equation (3) can be written into a simple analogy of the two-dimensional 
radial Schrodinger equation as [2,6] 



(E - V(r)) - w ^' ' \ R(r) = 0, (4) 

with the parameter 



V = t+\(D-2), (5) 

which is a linear combination of the spatial dimensions D and the angular mo- 
mentum quantum number £. 

2.1 The pseudoharmonic potential 

This potential has been studied in three dimensions [8] using the polynomial so- 
lution method. It has the following form [8] 



V(r) = D e (- 



r \ 2 



(6) 
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which can be simply rewritten in the form of isotropic harmonic oscillator plus 
inverse quadratic potential [9] as 

b 

V(r) = ar 2 + — + c, (7) 

where a = D e r~ 2 , b = D e r 2 e and c = —2D e . We want to solve the radial SE, Eq. 
(4), with the potential given in Eq. (7) by taking an ansatz for the radial portion of 
wave function 

R(r) = exp [p (a, r)] £ a n r 2n+s ^ 2 , (8) 

n=0 

where 

p(a,r) = -ar 2 . (9) 

Substituting Eq. (8) into Eq. (4) and equating the coefficient of r 2n+<5+3 / 2 to zero, 
one can obtain 

A n a n + B n+ ia n+ i + C ra+ 2a n +2 = 0, (10) 

where 



A n = -£(E-c) + 2a(2n + 5 + 2), (11) 



B„ = — 



2fi 

h 

2 fib 



h 2 



rj 2 -^ + (2n + S+^j (2n + 5 + ±y (12) 

C n = 0, (13) 



and the value of the parameter for p (a, r) can be evaluated as 



Q = ± Vir <14) 

To obtain a well-behaved solution at the origin and infinity, we must set a = 
~\f^~ wn i cn ensures that wave function ansatz, Eq. (8), be finite for all r. Fur- 
ther, for a given p, if a p ^ 0, but a p+i = a p+2 = a p+3 = • • • = 0, we then obtain 
A p = from Eq. (11), i.e., 



E s = c + 



2h 2 a 



(2p + 5 + 2). (15) 
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Carrying through a parallel analysis to Ref. [6], A n , B n and C n must satisfy the 
determinant relation for a nontrivial solution 



det 



B C\ ■ ■ ■ 

A B\ C2 








A v _ x B p 



0. 



(16) 



To utilize this method, we present the exact solution for p = 0, 1 as follows. (1): 
When p = 0, we can obtain from Eq. (15), the exact energy spectrum as 



E° = c + 



2h 2 a 



(6 + 2) 



(17) 



Further, it is shown from Eq. (16) that B = 0, which leads to the following 
restriction on the parameter 5 of the ansatz of wave function in Eq. (8) with r) and 
potential parameters as 



It follows that the wave function for p = can be written as 

^°\r)=a r s+2 - D / 2 exp 



1 2fm 2 

2V n 2 



(18) 



(19) 



where the normalization constant 



a = 



'2/ia 



8/2+1 



(5 + 1)! \ V h 2 

is determined from the requirement that 



(20) 



r°- x dr = 1. 



(21) 



(2): When p = 1, the exact energy spectrum becomes 

E{ = c + 



2h 2 a 



(5 + 4) 



(22) 



Furthermore, it is shown from Eq. (16) that B B 1 = A Ci, which leads to the 
restriction on 5 with the parameters of the potential and rj as 



(S + l) 2 -^- 



2/M ( ( ,,o\2 2 



(23) 



The wave function for p = 1 can be read as 



exp 



1 2p,a 2 



(24) 



where a and a\ are normalization constants. The relation between them can be 
determined by Eqs (10)-(14) as 



'§?<>,+ (<* + 3) 



which provides 



1 



2 2 /i6\ 
, T a 1 = 0, 



(5/2+2 



(25) 



'2/xa 



(5 + 2)V(<J + l)! IV ^ 



(26) 



Following this way, we can generate a class of exact solutions through setting 
p = 0, 1, 2, • • • , etc. Generally speaking, if a p ^ 0, a p+ \ = a p+2 = • • • = 0, we 
have A p = 0, from which we can obtain the energy spectra (cf. Eq. (15)). The 
wave function can be read 



^»( r ) = ( ao + a.r 2 + ■■■ + a p r 2p ) r 5+2 ~ D ^ 2 exp 



1 \2jm 2 
"2V ft 2 



(27) 



where (i = 0, 1, 2, • • • , p) are normalization constants. 

(i) For the pseudoharmonic potential [8], with the parameters following Eq.(7), 
the exact energy spectra become 



E n i = —2D e + A 



n 2 D, 



4n + 2 + 1 / 8 ^— - ^ 



Tr 



+ [2i+{D-2)\ 



,n,£ = 0,l,2,- 
(28) 



and radial wave function becomes 
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^ n \r) = (a + air 2 H h a„r 2n ) r' 



x exp 



1 


2fiD e 2 




u 2 r e 2 r 



(29) 



(ii) For the three dimensional (D = 3) isotropic harmonic oscillator plus inverse 



quadratic potential [9], a = hfj,u 2 , b = g and c = 0, the exact energy spectra are 



flUl 



n(' 



4n + 2 + J(2£+l) 2 + |^ 



,n,£ = 0,l,2,---,s>0 (30) 



and the corresponding radial wave function in this case reads 



^H(r) = (a + ai r 2 + ■■■ + a n r 2n ) r 



1+5 i+,/(2m) 



\2 i «M9 



cxp 



2a 



(3D 

where with i = 0, 1, 2, • • ■ , n are normalization constants and denotes the 
oscillator frequency. Its worthwhile to note that these results correspond to the 
results obtained by Ref. [9]. 

(iii) Also for a = B 2 and b = c = 0, the potential in Eq. (7) turns to the 
anharmonic oscillator potential. And then its exact energy spectra are given as [8] 

(32) 



E, 



h 2 



> —B (An + 2£ + D) , n, i = 0, 1, 2, 
\ 2/4 



and radial wave function becomes 



^ n \r) = (a + air 2 H h a n r 2n ) / exp 




Br 2 



(33) 



where a { with i = 0, 1, 2, • • • , n are normalization constants. On the other hand, 
setting a = /i = l,a = | and 6 = c = 0, gives the results of Ref. [22] for the 
results of exact harmonic oscillator energy states (cf. Eq. (8)) and wave functions 
(cf. after Eq.(4)). 
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2.2 The Mie-type potentials 

This potential has been studied in the D dimensions using the polynomial solution 
method [14]. An example on this type of potentials is the standard Morse [10] or 
Kratzer-Fues [11,12] potential of the form [10-12] 



' 2r P r\ " 
7i 



vir) = -De\=f (34) 



where D e is the dissociation energy between two atoms in a solid and r e is the 
equilibrium internuclear seperation. The standard Kratzer potential is modified 
by adding a D e term to the potential. A new type of this potential is the modified 
Kratzer- type of molecular potential [15] 



x 2 

r _ Te N 



V(r)=D e [— r Z) , (35) 

and hence it is shifted in amount of D e . The potential in Eq. (35) will be studied 
in D dimensions by making the wave function ansatz [6]. However, this potential 
has also been discussed before in three dimensions [15] and in D dimensions [14]. 
This potential [13] can be simply taken as 

a b 

V(r) = - + - + c , (36) 

where a = —D e r e , b = D e r\ and c = D e [15]. 

We take the following ansatz for the radial wave function 

R(r) = exp [p (a, r)] ]T a n r n+5+l '\ (37) 

n=0 

where 

p (a, r) = ocr. (38) 

Substituting Eq. (36) into Eq. (4) and setting the coefficient of r n+<5-1 / 2 to zero, 
we have 

A n a n + B n+ ia n+ i + C n+ 2d n +2 — 0, (39) 



where 



A n = -^ + 2a(n + 5+±y (40) 
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B 



c n = 0, 



and the value of the parameter for p (a, r) can be evaluated as 



« 2 = ~~ c )' 



from which we must set 



a — —\ 



2/i 

h 



2K E - c),E < c. 



(41) 
(42) 

(43) 
(44) 



which ensures that wave function ansatz Eq. (37) be finite for all r. This is required 
by the physically acceptable solution. On the other hand, for a given p, if a p ^ 0, 
but a v+ \ = a p+2 = a p+ 3 = • • • = 0, it is easy to obtain A p = from Eq. (40), i.e., 



Ei 



c — 



2% 2 [p + 5+\y 



(45) 



To apply this method, we will give the exact solutions for p = 0, 1 below. 

(1): When p — 0, it is found from Eq. (45) that the exact energy spectrum 
becomes 



En 



C — 



2ft 2 (*+ |) 



2 ' 



(46) 



The restriction on the parameter 5 with potential parameters and rj can be obtained 
from B = as 



2fib 

The wave function for p = now becomes 



(47) 



V> (0) (r) 



(P-2) 



a r 



cxp 



,E<c, 



(48) 



where 
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is a normalization constant obtained via Eq. (21). 

(2): When p — 1, the exact energy spectrum becomes 



s+i 



(49) 



E{ = c 



jj,a 



2K z (5+% 



2 ' 



(50) 



It is shown that the restriction on the parameter 5 with of the potential parameters 
and rj can be obtained from B B 1 = A C 1 



h 2 



The wave function for p = 1 becomes 



0. 



ip^\r) = (a + air) r 



5-(D-2)/2 



exp 



2/i 



2 (E - c)r 



,E <c, 



where the relation between them can be found as 

4/xa (8 + 1 



h 2 \8 + l 



a + (25 + 1) ai = 0, 



which provides 



(51) 



(52) 



(53) 



5+1 



Sj^a (f + 1) (J-%iE-6)\ .£■<,-. (54, 

R (2<S + 3) (26 + 1) v /(2i + 1)! V » R 

Following this way, we can generate a class of exact solutions through setting 
p = 0, 1, 2, • • • , etc. Generally speaking, if a p ^ 0, a p+ \ = a p+2 = • • • = 0, we 
have A p = 0, from which we can obtain the energy spectra (cf. Eq. (45)). The 
corresponding wave function can also be read as 



^(P\r) = (a + a x r + ■ ■ ■ + a p r p ) r 
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x exp 



2/i 



2 (E - c)r 



,E <c 



(55) 



where Oj (i — 0, 1, 2, • • • , p) are normalization constants. 

(i) For the standard Kratzer-Fues (Mie-type) potential [14], with the parame- 
ters a = —D e r e , b = D e r 2 and c = 0, the exact energy spectra are given as (cf. 
Ref. [14]) 



E - h 



h 2 



2n+ 1 



n 2 



+ [2i + D-2\ 



, (56) 



and the corresponding radial wave function becomes 



^(™)(r) = (ao + a±r H h a n r n ) r 



(D-2) 



_1 f±£e 
2 V 



erS +(2^+L>-2) 2 



x exp 



, i? n £ < 0. 



(57) 



(ii) For the modified Kratzer potential [15], with the parameters following Eq. 
(36), the exact energy spectra are given as (cf. Ref. [15]) 



E n l — D e — 



2// 



ft 2 



2n + 1 + 



a 2 



+ [2l + D-2\ 



(58) 



and the radial wave function becomes 



^(™)(r) = (a + air H h a n r n ) r" 



.(^ + i v /5^ + (2^Z)-2) 2 



x exp 



'-^ Ene -D e )r 



(59) 



(iii) Also after setting a = —A and b = c = 0, the potential in Eq. (36) turns to 
the Coulomb potential. Then its exact energy spectra are given as [23] 



E n e — 



M 2 



2h 2 (n + £+^) 



2 ■ 



(60) 
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and thus the radial wave function reads 



^H(r) = (a + air 2 + • • • + a n r") / exp 
where a { with i = 0, 1, 2, • ■ • , n are normalization constants. 

3 Concluding Remarks 

We have easily obtained the exact bound state solutions of the D dimensional 
radial Schrodinger equation for a diatomic molecule with two general potential 
forms representing the pseudoharmonic [8] and modified Kratzer molecular [15] 
potentials by using the wave function ansatz method [6]. The presented procedure 
in this study is systematical and efficient in finding the exact energy spectra and 
corresponding wave functions of the Schrodinger equation for various diatomic 
molecules. This new method is tested to calculate the energy spectra of the N 2 , 
CO, NO and CH molecules and our numerical calculations are similar to those 
given in Table II with parameter values shown in Table I (cf. Ref. [8]) for pseu- 
doharmonic potential and also to the ones given by Ref. [15] with the modified 
Kratzer potential for the given quantum numbers n = 0,1, 2, 3, 4, 5. This method 
is simple and promising in producing the exact bound state solution for further an- 
harmonic potentials [6], quarkonium potentials [21] and inter-nuclear potentials 
[24]. 
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